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Abstract 

The genome of the lactic acid bacterium Lactobacillus plantarum WCFS1 reveals the presence of a rich repertoire of 
esterases and lipases highlighting their important role in cellular metabolism. Among them is the carboxylesterase LpEstl a 
bacterial enzyme related to the mammalian hormone-sensitive lipase, which is known to play a central role in energy 
homeostasis. In this study, the crystal structure of LpEstl has been determined at 2.05 A resolution; it exhibits an afi- 
hydrolase fold, consisting of a central (3-sheet surrounded by a-helices, endowed with novel topological features. The 
structure reveals a dimeric assembly not comparable with any other enzyme from the bacterial hormone-sensitive lipase 
family, probably echoing the specific structural features of the participating subunits. Biophysical studies including 
analytical gel filtration and ultracentrifugation support the dimeric nature of LpEstl. Structural and mutational analyses of 
the substrate-binding pocket and active site together with biochemical studies provided insights for understanding the 
substrate profile of LpEstl and suggested for the first time the conserved Aspl 73, which is adjacent to the nucleophile, as a 
key element in the stabilization of the loop where the oxyanion hole resides. 

Citation: Alvarez Y, Esteban-Torres M, Cortes-Cabrera A, Gago F, Acebron I, et al. (2014) Esterase LpEstl from Lactobacillus plantarum: A Novel and Atypical 
Member of the «(? Hydrolase Superfamily of Enzymes. PLoS ONE 9(3): e92257. doi:10.1371/journal.pone.0092257 

Editor: Fernando Rodrigues-Lima, University Paris Diderot-Paris 7, France 

Received December 17, 2013; Accepted February 19, 2014; Published March 24, 2014 

Copyright: © 2014 Alvarez et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported by the Ministerio de Economi'a y Competitividad (BFU2010-17929/BMC) and the Factorfa de Cristalizacion (Consolider- 
lngenio-2007) (to Jose M. Mancheno), and AGL201 1-22745 and FUN-C-FOOD Consolider 25506 (to Rosario Munoz). The funders had no role in study design, data 
collection and analysis, decision to publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 

* E-mail: xjosemi@iqfr.csic.es 

9 These authors contributed equally to this work. 



Introduction 

Hydrolases constitute a class of enzymes that catalyse the 
hydrolysis of a wide variety of substrates, from peptides, amides or 
halides in addition to esters and triglycerides, as well as non- 
natural substrates. Although this assortment of substrates has 
complicated their classification, de facto they have been typified 
according to their known specificity. Esterases (EC 3.1.1), for 
instance, were defined as enzymes that hydrolyse ester bonds. 
Characteristically, they show specificity for either the alcohol or 
the acid moiety of the substrate, but not for both. Carboxyles- 
terases (EC 3.1.1.1), in particular, catalyses the hydrolysis of small 
carboxylic acid ester-containing molecules at least partially soluble 
in water, while lipases (EC 3.1.1.3) display maximal activity 
against water-insoluble long-chain triglycerides [1]. Thus, al- 
though catalytically similar, lipases and carboxylesterases must 
deal with physicochemically distinct environments: whereas lipases 
have to be capable of identifying an insoluble or heavily 
aggregated substrate, i.e. a water-substrate interface [2,3], 
carboxylesterase activity is maximal against monomeric substrates. 
Within this latter group of carboxylesterases substrate specificities 



vary widely, with some enzymes displaying highly specific activity 
towards particular esters such as acetylcholinesterase [4], heroin 
esterase [5] or Brefeldin A esterase [6], whereas others have 
activity against a broad range of substrates [7]. This group of 
enzymes is attractive for industry and in fact many carboxyles- 
terases have been utilized in the synthesis of ester compounds in 
non-aqueous solvents and also in stereospecific hydrolysis since 
they combine a broad specificity range with a high stereoselectivity 
[8-10]. 

The ESTHER database of lipases and esterases classifies these 
enzymes into four blocks, C, H, L and X [1 1] according to their 
amino acid sequence. Block H includes the hormone-sensitive 
lipase (HSL) family, a group of lipases and carboxylesterases from 
diverse biological sources which share sequence similarities with 
mammalian HSL [12]. Apart from the characteristic GXSXG 
motif around the active site serine, which is also found in serine 
proteases [13], they contain a highly conserved sequence of 
HGGG upstream the catalytic site. 

From a structural viewpoint, the amino acid sequence data on 
carboxylesterases indicate that they belong to the afihydrolase 
superfamily of enzymes [14—16]. Members of this superfamily 
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share a characteristic afl fold, which is based on an eight-stranded 
mostly parallel flsheet surrounded on both sides by a-helices. This 
structural framework supports a catalytic machinery based on a 
catalytic triad made up of a nucleophile, usually serine, an acid 
(Asp/Glu) and a histidine. The nucleophile is located within the 
above mentioned conserved G-X-S-X-G motif in a sharp turn 
between strand fS5and helix a3called the nucleophile elbow, where 
it can be approached by the substrate and by the hydrolytic water 
molecule [16]. The specific structural arrangement of the active 
site, particularly the geometry of the nucleophile elbow and the 
connecting loop in which the HGGG motif is located, contribute 
to the formation of the so-called oxyanion hole, usually formed by 
two or three backbone nitrogen atoms, which is needed to stabilize 
the negatively charged transition state [17]. 

Recently, crystal structures of several bacterial HSL enzymes 
have been reported, including carboxylesterases PestE and AFEST 
from the archaea Pyrobaculum calidifontis [18] and Archaeoglobus 
fulgidus [19], respectively, carboxylesterases from environmental 
samples EstE7 (PDB entry: 3k6k) EstE5 [20] and Est25 [21], 
esterases EstE2 and Sto-Est from Alicyclobacillus acidocaldarius [22] 
and Sulfolobus tokodaii [23], heroin esterase from Rhodococcus sp. [5], 
Brefeldin A esterase (BFAE) from Bacillus subtilis [6] and 
carboxylesterase Cest-2923 from Lactobacillus plantarum WCFS1 
[24]. A notable source of structural variability of this group of afi 
hydrolase enzymes emerges from the diverse quaternary structures 
observed: from monomers like EstE5 [20] to octamers identified 
within the crystals of AFEST [19]. These findings reveal the 
remarkable versatility of the ot(3 hydrolase fold in forming different 
assemblies, namely the presence of pleomorphism within this 
family of enzymes. An interesting example of this behaviour has 
been recently shown for Cest-2923, which may exist as mono- 
meric, dimeric and tetrameric species [24]. Although recent 
reports have focused on the analysis of quaternary assemblies 
within bacterial HSL enzymes [18,24], a systematic classification 
remains to be done. 

Here, we present the crystal structure of the carboxylesterase 
LpEstl from L. plantarum WCFS1 at 2.05 A resolution. This 
structure revealed that although LpEstl belongs to the bacterial 
HSL family of enzymes, as otherwise expected from its sequence, it 
represents a novel topological variant due to the presence of large 
and specific structural features that decorate a canonical afihy- 
drolase core. Additionally, the esterase forms dimers in the crystal 
and in solution through an association mode not observed yet in 
any other afShydrolase enzyme. Finally, in silico analyses of putative 
complexes revealed insights into the substrate specificity, which 
agree with the biochemical characterization of the enzyme. 

Results and Discussion 

Structure determination of LpEstl 

Crystals of recombinant, His-tagged LpEstl were prepared as 
previously described [25]. Analysis of the collected diffraction data 
revealed that these crystals were perfecdy twinned, belonging to 
the tetragonal space group 74, although they exhibited apparent 
point group symmetry 422 [25]. Unexpectedly, some of these 
crystals, when manipulated, spontaneously broke up into two 
parts, usually yielding two equal, untwinned crystals according to 
intensity statistics. This phenomenon not only provided a 
straightforward explanation for the identified perfect twinning 
and the apparent 422 point group symmetry (Figure SI), but also 
opened up the possibility to prepare untwinned crystals (or crystals 
with a very low twin fraction) suitable for structural studies. 
Optimized crystallization conditions found for His-tagged, seleno- 
methionine (Se-Met) labelled LpEstl were 1 M sodium malonate, 



0.5% (v/v) Jeffamine ED-2001, 100 mM HEPES, pH 7.0, and 
5 mM DTT (2 ul of protein 9 mg ml 1 plus 2 u.1 of reservoir 
solution). Despite the fact that both native and Se-Met labelled 
protein variants crystallized in different conditions, crystals of Se- 
Met LpEstl were essentially identical to the native ones, which 
allowed the collection of high quality, untwinned diffraction data 
(Table 1) that has permitted the determination of the structure of 
LpEstl by anomalous diffraction methods. The initial atomic 
model determined was then used as a molecular-replacement 
search model to obtain phases for a higher resolution data set 
(2.05 A resolution). This high-resolution model contains four 
independendy refined complete protein molecules (337 amino 
acids each) plus residues from the N-terminal TEV recognition 
sequence (4 in chains A, B and D, and 2 in chain C). Also, a total 
of 2196 water molecules, 2 malonate molecules and 5 glycerol 
molecules were modelled. The final refined model has an R w()r k of 
12.2% and an R free of 14.2%. The average B factor of the 
structure is 30.3 A . The analysis of Ramachandran plot showed 
that most of the modelled residues were in preferred and allowed 
regions (Table 1). 

LpEstl is a member of the bacterial HSL family 

The globular structure of LpEstl subunits has approximate 
dimensions of 45 Ax45 Ax50 A, and the overall secondary 
structure is a mixture of [S-sheets (20%) and a-helices (29%). 
The protein belongs to the afS-hydrolase superfamily [14] and is 
related to the bacterial hormone-sensitive lipase (HSL) family of 
enzymes [12] as concluded from structural similarity searches with 
DALI [26] . In this sense, multiple amino acid sequence alignment 
with CLUSTALW (Figure 1) reveals the presence in LpEstl of 
the sequence motif HGGG(A) that is known to contribute to the 
oxyanion hole [21]. As shown in Figure 1 this motif is localized 
within a sequence stretch highly conserved in the bacterial HSL 
enzymes, which encompasses part of strand [S5 and the connecting 
loop between this strand and helix oc3. The homologs with the 
highest structural similarity (Table 2) are hyperthermophilic 
carboxylesterases, namely the carboxylesterases PestE and AFEST 
from the archaea Pyrobaculum calidifontis [18] and Archaeoglobus 
fulgidus [19], respectively, followed by the mesophilic carboxyles- 
terases from environmental samples EstE7 (PDB entry: 3k6k) and 
EstE5 [20], and the thermophilic esterases EstE2 and Sto-Est from 
Alicyclobacillus acidocaldarius [22] and Sulfolobus tokodaii [23], respec- 
tively. Additionally, high structural similarity is also found with 
heroin esterase from Rhodococcus sp. [5] and Brefeldin A esterase 
(BFAE) from Bacillus subtilis [6]. Despite a low level of average 
sequence identity between LpEstl and its homologs (~24%), the 
3D structural similarity within this set of proteins is high (overall 
rmsd 2.3 A for ~280 Cot aligned atoms). This is mainly due to the 
presence of a common core [3-sheet surrounded by a-helices, 
which defines the canonical afi-hydrolase fold, which in turn 
typifies esterases with a Ser/Cys-His-Asp/Glu catalytic triad. On 
the contrary, high structural variability is observed in the so-called 
cap region situated over the active site on the carboxy-edge of the 
core fi-sheet (see below), specific for the bacterial HSL family of 
enzymes. 

Overall structure of LpEstl: a new variant of the ap 
hydrolase fold 

LpEstl consists of 337 residues with a molecular weight of 
-36.7 kDa (UniProtKB/Swiss-Prot code: Q88Y25). As a member 
of the afi-hydrolase superfamily, the structure exhibits a three- 
layered architecture, namely a central, almost parallel ten- 
stranded (5-sheet (two extra fistrands than the canonical fold) 
surrounded by five helices, two in the concave side of the sheet (a3 
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Table 1. Data-collection and refinement statistics. 







Native LpEstl 


SAD 


High resolution 


PDB code 


4c88 


4c87 


4c89 


Beamline 


ID14-4 (ESRF) 


ID29 (ESRF) 


ID14-4 (ESRF) 


Crystal parameters 


Space group 


14 


14 


14 


Unit-cell parameters (A) 


a = b 


168.34 


169.28 


168.79 


c 


184.20 


1 84.75 


184.57 


Matthews coefficient (A 3 Da~ 1 ) 


4.39 


4.44 


4.41 


Solvent content (%) 


72.0 


72.3 


72.1 


Data-collection statistics 


Wavelength (A) 


0.94000 


0.97915 


0.97914 


Resolution (A) 


53.68-2.65 


47.77-2.65 


53.82-2.05 




(2.79-2.65)" 


(2.79-2.65) 


(2.16-2.05) 


Total reflections 


615388 


1047518 


1207617 


Unique reflections 


74256 (10842) 


75273 (10961) 


161172 (23522) 


Emerge 


0.158 (0.508) 


0.111 (0.588) 


0.143 (0.702) 


Mean l/a 


12.1 (4.4) 


25.8 (4.9) 


10.1 (2.9) 


Completeness (%) 


100.0 (100.0) 


100.0 (100.0) 


100.0 (100.0) 


Anom. completeness {%) 




100.0 (100.0) 




Redundancy 


8.3 (8.3) 


13.9 (13.1) 


7.5 (7.5) 


Anom. redundancy (%) 




7.0 (6.6) 




Refinement statistics 


Protein molecules per ASU 


4 


4 


4 


Residues 


1352 


1356 


1365 


Waters 


1712 


823 


2196 


Glycerols 




17 


5 


Malonate 






2 


Total No. of atoms 


12192 


11502 


12839 


Rwork/Rfree (%) 


13.71/17.49 


14.91/18.25 


12.2/14.2 


Average fi factors (A 2 ) 


All atoms 


27.5 


23.8 


30.3 


Protein 


25.4 


22.9 


27.0 


Waters 


40.3 


29.7 


45.9 


Ligands 




66.6 


54.1 


R.m.s.d. from ideality 


Bonds (A) 


0.006 


0.007 


0.006 


Angles (°) 


0.960 


1.049 


1.025 


Ramachandran plot statistics 


Preferred regions (%) 


97.01 


97.80 


97.52 


Allowed regions (%) 


1.99 


2.12 


2.26 


Outliers 


0.00 


0.07 


0.22 



a Values in parenthesis are for the outermost resolution shell. 
doi:1 0.1 371 /journal.pone.0092257.t001 



and a8) and three in the convex side (oc4, 0(5 and oc7) (Figure 2A). 
The strand order for the core [3-sheet is (33, 04, 06, (35, |37, 08, 09, 
p 1 0, 02 and (51 with (51 and (34 being antiparallel to the others. 
The remarkable finding that this arrangement, where two N- 
terminal strands associate to the C-terminal, outermost strand of 
the core (3-sheet, has not been observed yet in any other 0(0- 
hydrolase enzyme confers a novelty to the molecular topology of 



the bacterial HSL family (Figure 2B). The overall quality of the 
2-Fobs - ^caic electron density map can be seen in Figure 2C. 

Together with this (3-sheet sandwiched by layers of helices, 
members of the bacterial HSL family show an additional, mainly 
helical cap domain [6,19,22]. The cap domain of LpEstl is formed 
by two sequence regions (residues 30-70 and 202-273), that group 
together around the carboxy-edge of the central |3-sheet. Here, 
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PI 



p2 



ctl 



ct2 



LpEstl 

PestE 

AFEST 

Sto-Est 

EstE7 

EstE5 

EstE2 

HEst 

BFAE 



1 10 20 30 40 50 60 

MPTINSIQTTVNGWKIVKPFNNDIAGEQFDPHVLQTLTAFKQPAILENDLAALRSGSLTPAIADPVGD 69 

MPLSPILRQILQQLAAQLQFRP-DMDVKTVREQFEKSSLILVK — MANEPIHRVEDI 54 

MLDMPIDPVYYQLAEYFDSLPKFDQ-FSSAREYREAINRIYEERNRQLSQHERVERVEDR 59 

MIDPKIKKLLESTIQLPIGKASVEEIRSLFKQFSSLTPR EEVGKIEDI 48 

MASMTGGNNMGRGSMGAMDQEIGTVTDTKMDPRDFLQLLKINAEKAEKNLPLDQKRAGMEALCERFPRA 65 

MASMTGGQQMGRGMAG PEIVKLK KILREKAVPPGTEVPLDVMRKGMEKVAFKAADD 43 

MPLDPVIQQVLDQLNRMPAPDYKHLSAQQFRSQQSLFPPVKKE PVAEVREF 52 

MTTFPTLDPELAAALTMLPKVDFADLPNARATYDALIGAMLADLS FDGVSLREL 54 

MTVTPERSTYTPPGRLGDESSGPRTDPRFSPAMVEALATFGLDEVAAAPPVSASDDLPTVLAAVGASHDGFQAVYDSIAL 80 



LpEstl 

PestE 

AFEST 

Sto-Est 

EstE7 

EstE5 

EstE2 

HEst 

BFAE 



P3 



(54 



P5 



a3 



P 6 , 



100 
HTvtVYFiTS 



80 90 

AVTVQSRNITALN RTVSVEWLTPQNVIN 

TIPGRGGP IRARVYRP-RDGERLPAyVYYl 

T I KGRNGD IRVRVYQQ 

TIPGSETN IKARVYYP 

EGVELTLTDL 



110 
FYGGVPGN- 



SFVLGSVET- 



-PVLVYY 



-KPDS- 

-KTQGPYGVtVYY^ 



SFVICSIES- 



SFVLGDIES- 



-GGVPC I RQAT DGAGAAH I L YF^ ^GY I S GSP ST 



IQVEQVTVAG CAAEWVRAPGCQAGKA- ILYL - 

DMDLPGRT LKVRMYRPEGVEPPYPALVYxJ 

SAPGLDGDPE VKIRFVTPDNTAGPVPVLLWI 



DLPTDRDDVETSTETILGVDGNEITLHVFRPAGVEGVLPGLVYT ■ 



120 130 

NjTVLLKLVAAKSHCEILNVDYS — 138 

HDHVCRRLANLSGAVWSVDYR — 116 

HDALCRRIARLSNSTWSVDYR — 119 

YDPLCRAITNSCQCVTISVDYR — 110 

HLVLTTQLAKOSSATLWSLDYR — 13 0 

5YVMGSINT — H RSMVGEISRASQAAALLLDYR — 107 

SWWGDLET — H DPVCRVLAKDGRAWFSVDYR — 115 

DPFCVEVARELGFAVANVEYR — 119 

161 



GFAIGTAES- 



3MTILTTDNRVHRRWCTDLAAAG-SVWMVDFRNA 



a4 



P7 



a5 



P8 



LpEstl 

PestE 

AFEST 

Sto-Est 

EstE7 

EstE5 

EstE2 

HEst 

BFAE 



140 1 
— LAPEAPAgAGILj 

— lapehkf^aave' 
— lapehrf^aavy 
— lapenkf^aaw 
— lapenpf^aavd 

— LAPEHPFgAAVE 
— LAPEHKFgAAVE' 
— LAPETTFgGPVN' 

wtaeghhpf"sgve^ 



50 160 170 

1laifqyleqrda etmitvagd 

\ydaakwvadnydklgvdngkiava5d 



180 190 200 210 

/IMAATNLNQQLGSNRINQQLLLYPVTAPNADHAG 

gGNLAAVTAIMARDRGESFVKYQVLIYP AVN 

i?GNLAAAVSIMARDSGEDFIKHQILIYP WN 

^GNLAAVTAILSK-KENIKLKYQVLIYP AVS 

cGGLTTASMLKAKEDGLPMPAGLVMLSP FVD 



]CYDATKWVAENAEELRIDPSKIFVGgD£ 
ALKWVYNNSEKFNGKYG-IAVGgD 

]C VAAYRALLKTAG SADRI I IAgD 

iVAAYRWLLDQGF KPQHLSISgL^.. ...... ....... 

alqwiaeraadfhldpariavgSdBaSgnlaavtsilakerggpalafqlliyp STG 

:yaallyihahaeelgidpsriavggqya5gglaagtvlkardegvvpvafqfleip eld 

claavlwvdehreslg — lsgvwq"e"g"gnlaiattllakrrgrldaidgvyasipyisggyawd 



213 
189 
192 
181 
189 
176 
188 
192 
241 



a7 



a8 



LpEstl 

PestE 

AFEST 

Sto-Est 

EstE7 

EstE5 

EstE2 

HEst 

BFAE 



220 



230 



240 



250 



260 



270 



280* 



290 



plwdlaafpiidsqrailtnyhdlfrqldsimtdyyvpenfdshsplisplhqenftmtgpttimvgefdpfrpqawayaqr 295 

ltgsptvsrveysgpeyviltadlmaw fgrqyfskpqdalspyas-pifadlsnli.palvitaeydplrdegelyahl 2 66 

fv-aptpsllef-geglwildqkimsw fseqyfsreedkfnplas-vifadlenlgpali itaeydplrdegevfgqm 267 

-fgqqylrsfadlldfrfs-piladlndlJJpaliitaehdplrdqgeayank 2 55 



FD- 



-LITKSLYDNGEGFFLTREHIDW- 



LT-LSRWSNSNLADRDFLAEPDTLGEM SELYVGG-EDRKNPLIS-PVYADLSGL^EMLIHVGSEEALLSDSTTLAER 27 7 

MT-CTNDSFKTRAEADPMVAPGGINKM AARYLNG-ADAKHPYAS-PNFANLKGLijPLLIHVGRDEVLLDDSIKLDAK 263 



YD-PAHPPASIEENAEGYLLTGGMMLW 



KLYAEA 265 



frdqylnsleelthpwfspvlypdlsgl^payiataqydplrdvg] 
dr-letvsmtnfvdtplwhrpnailswky-ylgesysgpedpdvsiyaapsratdltgij!ptylstmeldplrdegieyalr 272 
herrltelpslvendgyfienggmallvr-aydptgehaedpiawpyfa — sedelrglijpfwavneldplrdegiafarr 32 0 




300 310* 320 330 

LpEstl LAAADTATTFIQYQGLnHa — FA PLVD — QYWQSQDVAQVMAAALI 337 

PestE LKTRGVRAVAVRYNGVI^G — FVNFYPILE-EGREAVSQIAASIKSMAVA 313 

AFEST LRRAGVEASIVRYRGVL^G — FINYYPVLK-AARDAINQIAALLVFD 311 

Sto-Est LLQSGVQVTSVRFNNVI^G — FVSFFPFIE-QGRDAIGLIGYVLRKVFYGK 304 

EstE7 AGAAGVSVELKIWPDMPgV — FQMYGKFVN-AADISIKEICHWISARISKLAAALE 330 

EstE5 AKADGVKSTLEIWDDMIgV — WHAFHPMLP-EGKQAIVRVGEFMREQWAALAAALE 316 

EstE2 LNKAGVKVEIENFEDLI^G — FAQFYSLSP-GATKALVRIAEKLRDALA 311 

HEst llqagvsvelhsfpgtfkigsalvataavsergaaealtairrglrslspvs 323 

BFAE LARAGVDVAARVNIGLVi*5ADVI-FRHWLPAALESTVRDVAGFAADRARLRRP 372 



Figure 1. Protein sequence alignment between LpEstl and homologs from the hormone-sensitive lipase (HSL) family. PestE, 
carboxylesterase PestE from Pyrobaculum calidifontis; AFEST, carboxylesterase AFEST from Archaeoglobus fulgidus; Sto-Est, thermophilic esterase Sto- 
Est from Sulfolobus tokodaii; EstE7, esterase EstE7 from environmental samples; EstE5, esterase EstE5 from environmental samples; EstE2, thermophilic 
esterases EstE2 from Alicyclobacillus acidocaldarius; HEst, heroin esterase from Rhodococcus sp.; BFAE, Brefeldin A esterase from Bacillus subtilis. 
Residues from LpEstl in Residues from LpEstl in sterase fblue cylinders and orange arrows, respectively. Residues forming the catalytic triad are 
marked with an asterisk. Colour code for boxes is as follows: red, conserved residues in all proteins; yellow, highly conserved positions; cyan, residues 
that coincide with the expected canonical sequence motif characteristic of enzymes from the HSL family. 
doi:10.1371/journal.pone.0092257.g001 
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Table 2. Structural homologs of LpEstl as revealed by DALI. 





PDB ID 


Z score 


rmsd (A) 


Identity (%) 


NALI" 


NRES b 


PestE 


2YH2 


36.0 


2.1 


22 


289 


308 


AFEST 


1JJI 


34.6 


2.5 


26 


290 


311 


EstE7 


3K6K 


34.4 


2.4 


20 


285 


297 


EstE5 


3FAK 


33.9 


2.5 


18 


286 


297 


EstE2 


1EVQ 


33.3 


2.4 


25 


288 


308 


Sto-Est 


3AIK 


32.4 


2.2 


23 


269 


283 


Heroin esterase 


1 LZL 


32.4 


2.4 


26 


285 


317 


Brefeldin A esterase 


1 JKM 


32.0 


2.4 


21 


291 


361 



a NALI: number of aligned residues; 
b NRES: total number of residues. 
doi:1 0.1 371 /journal.pone.0092257.t002 



three oc-helices (ocl, oc2, and a7) and irregular loops are identified. 
Helix 0(7, situated between strands (38 and f39, is exceptionally long 
when compared to the equivalent ones from the rest of the 
bacterial HSL members (see below) and mainly interacts with the 
extended al-a2 connecting loop, which lies almost parallel along 
its axis. This helix, together with helices ocl and 0(2, form the upper 
walls of a funnel-like structure, which has the catalytic machinery 
at its base, with a surface area of 683 A as determined with the 
CASTp server [27]. Conversely, the lower walls of this funnel are 
contributed by the loops between strand f35 and helix 0(3 and 
between helix 0(8 and strand filO. Among the residues facing the 
funnel ~70% are hydrophobic, indicating the predominance of 
hydrophobic interactions in this environment. 



Structural comparisons of the cap regions of LpEstl and its 
bacterial HSL homologs reveal high structural variability 
(Figure 3). Nonetheless, despite this variability three distinct 
patterns can be easily identified, which are represented by LpEstl 
(Figure 3A) BFAE (Figure 3B) and the rest of the HSL enzymes 
(Figure 3C), respectively. This latter, predominant pattern reveals 
that the main source of structural variability resides in the different 
lengths and relative orientations of N-terminal helices (equivalent 
to helices al and 0(2 from LpEstl) and the corresponding 
connecting loop. In contrast, the C-terminal part of the cap 
region of these enzymes (equivalent to residues 200-273 from 
LpEstl) remains highly conserved (overall rmsd 1.73 A), even in 
BFAE whose the N-terminal part clearly departures from the 




Figure 2. Crystal structure of the LpEstl subunit. (A) Ribbon representation of the LpEstl subunit; two different views are depicted. Canonical 
P-strands forming the core (3-sheet are shown in yellow, whereas the two N-terminal, non-canonical ones are shown in orange. a-Helices are in blue, 
except the first two previous to the first p-strand from the core p-sheet (P3), which are shown in green. (B) Topology diagram of the LpEstl fold. 
Colour code for the secondary structure elements are as in (A). The positions of the residues forming the catalytic triad Ser174, Asp283 and His313 are 
indicated as S, D and H, respectively, and those forming the oxyanion hole, Glyl 07 and Ala108, are indicated as GA. (C) Representative 2F obs - F cak 
density map contoured at 1a. 
doi:1 0.1 371 /journal.pone.0092257.g002 



PLOS ONE | www.plosone.org 



5 



March 2014 | Volume 9 | Issue 3 | e92257 



Structure and Dimeric Assembly of Esterase LpEstl 




Figure 3. Three-dimensional comparisons of cap regions from bacterial HSL enzymes. (A) Superposition of the cap regions of LpEstl 
(yellow) and PestE from P. calidifontis (green), which has been chosen arbitrarily as a representative model of the predominant pattern (see the text). 
(B) Superposition of the cap region from BFAE (dark green) and that from PestE. (C) Superposition of cap regions from PestE (green), AFEST (cyan), 
EstE2 (grey), EstE5 (light brown), EstE7 (magenta), heroin esterase (pale green) and Sto-Est from (blue). All structures are represented as ribbon models. 
doi:1 0.1 371 /journal.pone.0092257.g003 



predominant pattern due to the relative displacement of helix oc2 
and the much longer otl— oc2 connecting loop. Conversely, the cap 
region of LpEst 1 differs from the main pattern both in the N- and 
C-terminal parts due to the different relative orientation of helices 
aland a2 and the much longer helix a 7 (22 residues versus an 
average of 10 residues for the rest of the proteins). 

Dimer arrangement 

The four independent molecules of the LpEst 1 asymmetric unit 
are arranged in two identical dimers, consisting of the pairs A— C 
and B-D. The independent molecules are almost perfectly 
superimposable with an overall rmsd value of 0.42 A as 
determined with PDBeFOLD [28]. According to the analyses 
with PISA [29] and PIC [30] servers, the subunits within each 
dimer generate an interface area of 1130 A and interact via 19 
hydrogen bonds and 24 hydrophobic interactions. The contacting 
interface is basically made up of amino acid side chains from the 
long connecting loops between helix a6 and strandfS9 and between 
helix a5 and strandf38, and also from residues contributed by 
helices 0t4and a5 from each subunit (Figure 4A). It is remarkable 
that this region of association between LpEstl subunits, and 
therefore its association mode, is different from those observed for 
the dimeric assemblies of the members of the HSL family, which 
involves the antiparallel association of the C-terminal, outermost 
fistrand from the core [3-sheet (herein referred as to the canonical 
strand [38) (Figures 4B and 4C), as well as interactions between 
additional structural elements, as indicated below. Indeed, this is 
the case for the LpEstl homologs, which are all (at least) dimeric 
species with the exception of the monomeric carboxylesterases 
EstE2 [22] and EstE5 [20]. Dimers adopting the above described 
general topological arrangement are herein defined as canonical 
ones. A detailed structural analysis of these canonical dimers 
reveals an underlying complexity, partially advanced before [18], 
that points to the existence of at least two distinct subtypes of 
dimers: subtype 1 canonical dimers, formed by PestE, AFEST, Sto- 
Est, and Brefeldin A esterase, and subtype 2 canonical dimers, 
formed by Est7 and heroin esterase. In subtype 1 dimers, together 
with the antiparallel association of the canonical strands (38 
(mainly through their C-terminal ends), a large dimerization 
interface is made up of amino acid side chains contributed by the 
C-terminal ot-helix from each subunit. In particular, interactions 
are primarily observed between the two participant C-terminal a- 
helices (mainly between the central part and C-terminal end of 
each helix) and also between the N-terminal end of these same a- 
helices and the N-terminal end of the canonical strand [38 from the 
other subunit (Figure 4B). Conversely, in subtype 2 canonical 
dimers the main dimerization interface is formed by amino acid 
side chains from the canonical strands [38 and the immediately 
upstream a-helices. In this case, the interactions are mainly 



established between equivalent structural elements from each 
subunit, that is, between the C-terminal, canonical (3-strands 
themselves and the participant a-helices (Figure 4C). A 
topological corollary that can be inferred from this structural 
classification is that each dimer subtype is unambiguously defined 
by the side of the subunits (with respect to the plane of the central 
[3-sheet) that participates in dimer formation together with the 
antiparallel association of the strands (38. If these sides were 
arbitrarily defined as cis and trans, respectively (being the cis side 
the one in which the a-helix after the canonical strand [38 is 
located), subtype 1 dimers can be defined as resulting from a cis-cis 
subunit association, whereas subtype 2 dimers would result from 
trans-trans subunit associations. It is notable that this structural 
classification correlates with the fact that proteins belonging to the 
subtype 1 of dimers form stable, higher order assemblies within the 
crystals, with the dimers as basic building blocks, whereas proteins 
from the subtype 2 do not. Thus, tetramers are observed in the 
four subtype 1 proteins considered (the biological assembly of 
Brefeldin A esterase has been assigned by the authors to be a 
tetramer; see PDB entry ljkm), and even octamers in the case of 
AFEST. Structurally, the observed tetramers result from the 
association of subunits through their trans sides (Figure S2). This 
suggests that both cis and trans sides may be involved in dimer 
formation but only trans sides would be involved in tetramer 
formation. Of course, this does not necessarily exclude the 
possibility of cis-cis tetramers but the current structural information 
suggests that, should they exist, they are not abundant or 
canonical. 

The novelty of the association mode for LpEstl subunits 
observed within the crystal prompted us to characterize the 
oligomeric state of LpEstl in solution. Thus, as a first approach, 
we estimated the molecular weight of the enzyme by chromato- 
graphic and ultracentrifugation analyses (Figure 5). Results from 
analytical gel filtration experiments revealed that LpEstl behaves 
in solution (20 mil Tris, pH 8.0, 0.1 M NaCl and 5 mM DTT) as 
a unique species with an apparent molecular weight of 78 ±3 kDa 
(n = 3) (Figure 5A), which compares well with the value 
theoretically expected for the dimer (77 kDa). Secondly, analytical 
ultracentrifugation studies revealed that LpEstl behaves in 
solution as a single, homogeneous species with a sedimentation 
coefficient of 4.5 S (sedimentation velocity experiments) whose 
behaviour, as derived from sedimentation equilibrium experi- 
ments, fitted well to an ideal model of a unique species with a 
molecular mass of 77.2±4.2 kDa (n= 3) (Figure 5B). Therefore, 
these studies support the dimeric character of LpEst 1 in solution. 

In silico and mutational analyses of the dimer interface 

Once the dimeric nature of LpEst 1 was demonstrated both in 
the crystal and in solution, we further analysed the contacting 
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A 




Figure 4. Dimeric assembly of LpEstl and classification of the dimers of the HSL family members. (A) Two orthogonal views of a dimer 
of LpEstl . Each subunit is shown as ribbon model with different colour. (B) Subtype 1 of dimers (PDB entry: 3aik), characterized by a cis-cis association 
of subunits. C/s is arbitrarily defined as the side of the protein with respect to the plane of the core (3-sheet where the a-helix downstream the 
canonical strand (38 is situated (see the text for details). (C) Subtype 2 of dimers (PDB entry: 1 1zl), characterized by a trans-trans association of subunits. 
Trans is arbitrarily defined as the side of the protein with respect to the plane of the core (3-sheet where the oc-helix upstream the canonical strand (38 
is situated. The orientation of is this dimer is as in (B). 
doi:1 0.1 371 /journal.pone.0092257.g004 



regions in terms of interaction energies for dimer formation of 
both LpEst 1 and, for comparison purposes, also of its canonical 
homologs (Figure 6). Some important conclusions can be derived 
from this analysis: first, LpEstl displays a low dimer stabilization 
energy relative to the rest of the dimeric proteins with most of this 
stabilization originating from apolar contacts, together with some 
highly directional hydrogen bonds; second, the carboxylesterases 
PestE [18] and AFEST [19] exhibit the highest stabilization, in 
agreement with their hyperthermophilic character, with the most 
important contributions arising from coulombic interactions and 
hydrogen bonds; third, a very high stabilization is observed for 
Brefeldin A esterase and heroin esterase, despite their mesophilic 
character; fourth, further stabilization is attained upon tetramer 
formation of PestE, AFEST and Sto-Est, with the main driving 
force for this association being hydrophobic interactions 
(Figure 6). The high stabilization for Brefeldin A esterase is 
probably explained by the fact that, in addition to the canonical 
contacting interfaces of subtype 1 dimers, there are numerous 
intersubunit contacts involving the large and protein-specific cap 
region. The explanation for the high stability of heroin esterase 
remains an open question. 

Regarding the contacting region of LpEstl subunits, energy 
decomposition pinpointed six residues whose calculated contribu- 
tions are larger than —6 kcal mol ': Leu81, Gin 189, Leu 190, 
Asn253, Phe254 and Leu260. The fact that these residues are 
distributed all along the contacting interface (Figure 7) suggests a 
homogeneous stabilization in this region. Almost in the geomet- 
rical centre of the interface, a bidentate hydrogen bonding 
interaction is identified between the side-chain carboxamide 
groups of Gin 1 89 from each subunit. Taken these characteristics 
into account, we raised the working hypothesis that incorporation 
of a charged residue at this position might destabilize this 
environment. With this aim, we have produced and purified the 
Glnl89Glu mutant. We observed that after purification (IMAC on 



a His-Trap FF column plus size exclusion chromatography) the 
protein resulted metastable in solution since it slowly aggregated if 
maintained in 20 mM Tris-HCl, pH 8.0, 0. 1 M NaCl, in contrast 
to the wild-type protein, which remained stable. However, when 
the last chromatographic step is developed under acidic conditions 
in Mcllvaine buffer pH 5.0 (Na 2 P0 4 , citric acid, pH 5.0) no 
aggregation was observed for the mutant variant when stored at 4 
°C for 2-3 days, which is dimeric as revealed by sedimentation 
equilibrium experiments in these acidic conditions (Figure S3). 
These results indicate the formation of metastable monomeric 
species under neutral conditions, which eventually aggregate upon 
the exposure of the hydrophobic interface to the bulk solvent. 
However, dimer formation is observed under acidic conditions due 
to the protonation of the Glu 1 89 side chain. Moreover, in these 
acidic conditions, the Glnl89Glu mutant exhibits a hydrolytic 
activity against /)-nitrophenyl acetate almost equivalent to that of 
the wild-type enzyme. 

In summary, taken as a whole, these results demonstrate that 
LpEstl is a dimer despite the relative low overall interaction 
energies calculated at the interface and further suggest, that the 
apolar interactions are the main driving force for dimerization. 

Active site and substrate binding pocket 

The structural homology found with enzymes from the bacterial 
HSL family permitted the straightforward identification of the 
catalytic machinery of LpEstl as a classical catalytic triad with 
Serl74 as the nucleophile, His313 as the general base that 
deprotonates the serine hydroxyl and Asp283 as the residue that 
increases the pK a of the histidine imidazole ring. We validated this 
result experimentally since the Serl74Ala, His313Ala and 
Asp283Ala mutants displayed no catalytic activity in terms of 
hydrolytic activity against /j-nitrophenyl acetate. Significant 
conformational changes affecting to the global fold of the proteins 
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Figure 5. Analysis of the oligomeric state of LpEstl in solution. 

(A) Analytical gel-filtration of LpEstl on Superdex 200 1 0/300 GL Tricorn 
column. The elution profile of LpEstl is shown together with the elution 
positions for some standard proteins (molecular mass in kDa). Inset, 
semilog plot of the molecular mass of all standards used versus their K av 
values (open circles). The closed circle indicates the position of the K av 
value of LpEstl interpolated in the regression line (solid line) (B) 
Analytical ultracentrifugation analysis of LpEstl. Sedimentation equilib- 
rium analysis of LpEstl (10 uM) in Tris buffer (20 mM Tris-HCI, pH 8.0, 
and 0.1 M NaCI) at 9,000 rpm (open squares) and 13,000 (open circles). 
Absorbance at 280 nm is plotted against the radial position from the 
center of the rotor. The fit to the data set (solid line curves) corresponds 
to an ideal species with a molecular mass of 77.2±4.2 kDa (n = 3). 
Residuals from this fit are shown in the panel at the bottom. 
Calculations were done with the program Heteroanalysis [47]. Inset, 
sedimentation coefficient c(s) distributions for LpEstl (10 uM) in Tris 
buffer (20 mM Tris-HCI, pH 8.0 with 0.1 M NaCI). Raw sedimentation 
velocity profiles for this analysis were acquired using absorbance at 
280 nm, 45,000 rpm, 20 °C, and different times (not shown). 
Calculations were done with the program Sedfit [46]. 
doi:1 0.1 371 /journal.pone.0092257.g005 

as a result of the included mutations could be discarded by far-UV 
CD measurements (not shown). 

The three catalytic residues are located at canonical positions 
within the Dtp hydrolase fold [14—16], at the carboxy-edge of the 
core P-sheet. The nucleophile Serl74 is found at the apex of the 
so-called "nucleophilic elbow" between strand (57 and helix 0(5, 
within the conserved sequence motif GX!SX 2 G (GDSAG), a 
signature of the hydrolase family [14—16]. The constrained 
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Figure 6. Intersubunit binding energy decomposition for dimer 
and tetramer formation. Colour code is as follows: black, electro- 
static term; orange, van der Waals term; blue, hydrogen bond term; 
yellow, total contribution, including desolvation energies (not shown). 
doi:10.1371/journal.pone.0092257.g006 

conformation of the nucleophile peptide backbone (tp = 56° and 
v[i = — 1 28°) facilitates the formation of the "oxyanion hole" 
responsible for the stabilization of the negatively charged 
tetrahedral intermediates of the catalytic reaction [16]. Structural 
comparisons with other bacterial HSL homologues suggest that 
the oxyanion hole of LpEstl would be formed by the backbone 
nitrogen atoms of Glyl07, Ala 108 and Ala 1 75 namely it would be 
a tridentate structure similarly to those of heroin esterase [5], 
esterase EstEl [31] or carboxylesterase EstE2 from Alkyclobacillus 
acidocaldarius [22]. Nonetheless, as will be shown below, the in silico 
analysis of this region discarded any role for Ala 175 in the 
oxyanion hole. 

A complex network of hydrogen bonds is identified in the active 
site (Figure 8). Thus, the Oyatom of Serl74 establishes a 
hydrogen bond (2.8 A distance) with the Ns2 atom of the 
imidazole ring of His313, situated in the 13-residue loop between 
strand (510 and helix 0(9. In turn, the N51 atom of this latter 
residue is at hydrogen bond distance to 051 (3.1 A) and 082 
(2.8 A) atoms of Asp283, this residue being further stabilized by a 
2.6 A hydrogen bond with a water molecule, prototypical for the 
members of the HSL family [12]. This solvent molecule also forms 
hydrogen bonds to the carbonyl oxygen of Asn312 (2.8 A) and the 
amide nitrogen of Gly280 (2.7 A). 

The detailed structural characterization of the LpEstl active site 
permitted us to carry out docking experiments with the substrates 
phenyl acetate, triacetin and tributyrin (see below Biochemical 
characterization) (Figure 9). The best poses obtained from the 
automated docking protocol were consistent with the expected 
orientation of an ester-containing substrate in an esterase active 
site: the carbonyl moiety of the ester group is located near the 
catalytic hydroxyl group from Serl74 and is basically stabilized by 
two hydrogen bonds from the NH groups of Glyl07 and Alal08, 
which would then form the oxyanion hole. This result rules out the 
backbone nitrogen atom of Ala 175 as a participant in the 
oxyanion hole, which would thus be bidentate [16]. In addition, 
these three enzyme:substrate complexes reveal that the acid 
moiety of the substrates is inserted into a small, hydrophobic 
subpocket (S 1) lined by Met245, Phe285, Leu24 1 and Val204, and 
the alcohol part would lie in a open and large subpocket (S2). 
Apparendy, this particular architecture of the substrate binding- 
pocket in two distinct subpockets points to the hydrolytic activity of 
LpEstl being directed against esters with small acid moieties but 
larger alcohol moieties. The analysis of substrate specificity 
correlates well with this prediction (see below). 
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Figure 7. Distribution of the energetically relevant residues within the LpEstl contacting interface. These residues are basically 
hydrophobic in agreement with the relevance of this type of interactions in the stabilization of the dimer. The bidentate hydrogen bonding 
interaction between the side-chain carboxamide groups of Gin 189 (close up view) is situated at the core of the interface. The 2F obs - F cak density map 
is contoured at 1 a. 

doi:1 0.1 371 /journal.pone.0092257.g007 



A notable feature of the surroundings of the nucleophile Serl74 
is the presence of the acidic residue Asp 173. This is unexpected 
considering the hydrophobic character of the side chains that 
make up the binding pocket in LpEstl. Remarkably, all its 
structural homologs also have an equivalent acidic residue except 
for heroin esterase (Figure 1). In LpEstl, Aspl73 participates in a 
network of polar interactions, some of them mediated by water 
molecules, and this makes the carboxylate moiety to be highly 
oriented towards the loop making up the oxyanion hole 
(Figure 10). Thus, the oxygen atom OSl is at hydrogen bonding 
distance to the amide nitrogen of Glyl06 (3.2 A), the hydroxyl 
group of Tyrl03 (2.7 A) and a highly ordered water molecule 
(3.0 A). Interestingly, the interaction between Aspl73 and Tyrl03 
is conserved in all HSL proteins considered herein (with the 
exception, again, of heroin esterase). Conversely, the oxygen atom 
082 is at hydrogen bonding distance from the same amide 
nitrogen of Glyl06 (3.1 A) and a solvent water molecule (3.0 A), 
which in turn interacts with the hydroxyl group of Tyr202 (2.7 A). 
The orientation towards the loop making up the oxyanion hole 
and the well-ordered character of the Asp 173 side chain suggest 
that this residue may play an important structural role in LpEst 1 , 
and presumably in this group of enzymes, particularly in the 
stabilization of the proper conformation of the oxyanion hole. In 




agreement with this hypothesis, the replacement of this residue by 
Ala (Aspl 73Ala) resulted in a folded (Figure S4) but fully inactive 
protein variant unable to hydrolyze /;-nitrophenyl acetate. It is 
obvious that the possibility that Asp 173 plays a direct functional 
role cannot be discarded since its acidic character should 
contribute significandy to the negative electrostatic potential 
around the active site, a feature displayed by lipases and esterases 
in the pH range associated with their maximum activity [2] . This 
aspect is currently under investigation. 

Biochemical characterization 

We have examined some important enzymatic properties of 
LpEstl (Figure 11). The optimum pH for hydrolytic activity 
against /;-nitrophenyl acetate is 6.5 (Figure 11A), which is a value 
typically observed for esterases, in contrast to the higher pH values 
(~8.0) displayed by lipases [2]. Regarding to temperature, the 
protein presented highest activity at ~30 °C, although at 37 °C 
exhibited a high level of activity (~50%) (Figure 11B). These 
values for optimum pH and temperature are commonly found in 
other esterases from Lactobacilli [32-34]. On the other hand, 
temperature stability measurements show a drastic reduction in 
LpEstl hydrolytic activity upon incubation of the esterase at 37 °C 
(Figure 11C) 




Figure 8. Stereoview of the network of interactions present around the catalytic triad of LpEstl. Residues forming the catalytic triad 
(Ser174, His313 and Asp283) and those coordinating a well-ordered water molecule that is also identified in the HSL family members {blue sphere) are 
shown as sticks. The 2F 0 t, s - F ca i c density map is contoured at 1a (green: amino acid side chains; blue: water molecule) Distances are in A. 
doi:1 0.1 371 /journal.pone.0092257.g008 
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Figure 9. Stereoviews of three putative LpEstl :substrate complexes resulting from docking studies with CRDOCK. The three 
complexes correspond to phenyl acetate (a), triacetin (b) and tributyrin (c). In all cases, interactions are observed between the nucleophile (Ser174) 
and also reveal the stabilizing effect of the oxyanion hole formed by the backbone nitrogen atoms of Gly107 and Ala108. Residues are shown as stick 
models. Distances are in A. 
doi:1 0.1 371 /journal.pone.0092257.g009 



The acyl-length selectivity against />-nitrophenyl ester substrates 
follows this order: C2>C4>C8>C12>C14, indicating a prefer- 
ence for short acyl-length esters (Figure 11D). The kinetic 
parameters for C2 and C4 substrates were determined spectro- 
photometrically. In both cases, LpEstl exhibited a hyperbolic 
Michaelis-Menten kinetics (not shown). The kinetic parameters are 



shown in Table 3. From the values of these parameters it can be 
deduced that the catalytic efficiency (k^,/ K^) for pNPA hydrolysis 
is around 8-fold the one observed for pNPB hydrolysis. 
Conversely, the study of the substrate profile has been qualitatively 
analysed with the use of a library of esters as described previously 
[35]. This study reveals maximum hydrolysis against phenyl 
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Figure 1 0. Stereoview of the environment around the residue Aspl 73. The side chain of the Asp1 73 residue is highly oriented towards the 
loop where the oxyanion hole resides probably contributing to its stabilization. The network of hydrogen bonds is shown as dashed lines. 
Participating residues are shown as sticks and water molecules as blue spheres. Distances are in A. 
doi:10.1371/journal.pone.0092257.g010 



acetate, and significant activity against triacetin, tributyrin and 
isopropenyl acetate (Figure S5). Apparendy, this activity profile 
for LpEstl basically suggests a preference for a small acid moiety 
of the substrate, e.g. acetate or butyrate, in good agreement with 
the close configuration of subpocket SI. Conversely, the absence of 
hydrolytic activity against other esters such as ethyl acetate, vinyl 
acetate and propyl acetate, for instance, clearly indicates an 
important influence on hydrolytic activity of the alcohol moiety of 



the ester. In this regard, the combination of the relaxed, open 
conformation of subpocket S2 with the absence of hydrolytic 
activity against esters with small alcohol moieties point to a 
threshold in the size (or length) of this part of the substrate that 
would be necessary to accomplish for the formation of a 
productive enzymeisubstrate complex. Hence, as a whole, these 
two aspects, namely narrow substrate profile and recognition of 
both the acid and alcohol moieties of the substrate, would qualify 




time (h) 



Figure 11. Biochemical characterization of LpEstl. (A) Dependence on pH of hydrolytic activity of LpEstl against pNPA. (B) Dependence on 
temperature of hydrolytic activity of LpEstl against pNPA. The optimum temperature for esterase activity was —30 r C. (C) Analysis of the temperature 
stability of LpEstl. The enzyme was incubated in 50 mM sodium phosphate buffer pH 7.0 at 22 °C (closed triangles), 30 °C (open triangles), 37 °C 
(closed circles), 45 "C (open circles) and 55 °C (closed squares) for 15 min, 30 min, and 1, 2, 4, 6 and 20 h. The values correspond to the mean of three 
independent experiments. (D) Dependence of the esterase activity of LpEstl on the chain length of p-nitrophenyl (p-NP): p-NP acetate (C2), p-NP 
butyrate (C4); p-NP caprylate (C8); p-NP laureate (CI 2); and p-NP myristate (C14). 
doi:10.1371/journal.pone.0092257.g011 
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Table 3. Kinetic parameters for pNPA and pNPB hydrolysis by LpEstl. 





Substrate 


V max (fimol min 1 mg 


K m [mM] 


kc (s"') 




mM _1 ) 


pNPA 


75±4 


0.38±0.05 


40.4±2.0 


1 06 ± 1 2 




pNPB 


8+1 


0.2+0.1 


4.4±0.3 


14.7±3.3 





Enzyme activities were determined at 30 °C in 50 mM sodium phosphate buffer, pH 7.0. Results are the mean value ± SD from three independent experiments. 
doi:1 0.1 371 /joumal.pone.0092257.t003 



LpEstl as a specific esterase, probably reflecting a distinct 
biological role of the enzyme, which in turn can be linked to its 
unique structural peculiarities. 

Materials and Methods 

DNA manipulations 

The expression vector pURI3-TEV- LpEstl coding for the wild- 
type esterase LpEstl from L. plantarum was used as template for the 
preparation of single point mutants affecting to the catalytic 
residues, Aspl73 and Glnl89 residues, essentially as previously 
described [25,36]: Serl74Ala, Asp283Ala and His313Ala, As- 
pl73Ala and Glnl89Glu. The mutagenic primers used (forward 
and reverse, respectively) were: (5 ' -C GTTGCTGGCGATGC - 
GGCTGGCGC-3') and (5'-GCGCCAGCCGCATCGCCAG- 
CAACG-3') for the Serl74Ala mutant, (5 ' -GGCGAGTTTGC - 
TCCCTTCCG-3') and (5'-CGGAAGGGAGCAAACTCGCC- 
3') for the Asp283Ala mutant, (5'-GCTTGAACGCCGCTT- 
TCGCACC-3') and (5'-GGTGCGAAAGCGGCGTTCAAGC- 
3') for the His313Ala mutant, (5 ' -GATTACC GTTGCTGGC- 
GCTTCGGCTGGCGC-3 ') and (5'-GCGCCAGCCGAAGCG- 
C C AGC AAC GGT AATC - 3 ' ) for the Aspl73Ala mutant, and (5'- 
GCAACAAATCTGAATCAAGAACTCGGTAGC-3 ') and (5'- 
GCTACCGAGTTCTTGAT TCAGATTTGTTGC-3 ') for the 
Glnl89Glu mutant. Production and purification of the proteins 
were as previously described [25] . 

Overexpression and protein production 

Expression and purification of Se-Met labelled LpEstl esterase 
was done essentially as described for the native protein [25], using 
a minimum medium containing selenomethionine and buffers 
supplemented with 5 mM DTT. Final production yield was 6 mg 
per litre of culture. Incorporation of Se-Met was checked by mass 
spectrometry on a Finnigan LCQDeca ion trap Mass Spectrom- 
eter (Thermo Electron, San Jose, CA, USA). 

Crystallization and data collection 

Diffracting, perfectly twinned crystals of native LpEst 1 esterase 
were prepared as previously described [25]. For data collection, 
20% (v/v) glycerol was added to the crystallizing precipitant as a 
cryoprotectant and the crystals were flash-cooled in a 100 K 
nitrogen gas stream. Unexpectedly, untwinned crystals could be 
prepared by mechanically manipulating the first ones (see below 
for details). A diffraction data set was collected on ID 14-4 
beamline at the ESRF (Grenoble, France) using a wavelength of 
0.9400 A. An ADSC Quantum Q315r CCD detector was used 
with a crystal-to-detector distance of 304.72 mm, and a total of 
200 images were collected with a 1° oscillation angle. The 
diffraction data were processed with iMOSFLM [37]. These 
untwinned crystals belong to the tetragonal space group 74 and 
have unit cell parameters a = b= 168.34 A, c= 184.20 A 
(Table 1). Conversely, optimized crystals of the Se-Met labelled 
esterase were prepared in a condition consisting of 1 M sodium 



malonate, 0.5% (v/v) Jeffamine ED-2001, 100 mM HEPES, 
pH 7.0, and 5 mM DTT (2 \ll of protein 9 mg ml -1 plus 1 |il 
of reservoir solution). Similarly to the above crystals prepared with 
unlabelled protein, these new crystals were perfectly twinned, and 
their manipulation also permitted the production of untwinned 
crystals (or twinned crystals with a low twin fraction). The SAD 
data set was collected on beamline ID29 at the ESRF (Grenoble, 
France) using a wavelength of 0.97915 A. The detector was a 
Pilatus 6 M. A total of 3600 images were collected with a 0.1° 
oscillation angle and a crystal-to-detector distance of 388.02 mm. 
The diffraction data were processed with the XDS [38] program 
package. The crystal belonged to the tetragonal space group 74, 
with unit cell parameters a = b= 169.28 A, c= 184.75 A 
(Table 1). Interestingly, these latter crystals on average diffracted 
at higher resolution than those prepared with unlabelled protein. 
Hence, the high-resolution data set was collected with one of these 
crystals. This data set was collected on beamline ID 14-4 at the 
ESRF using a wavelength of 0.97914 A and an ADSC Quantum 
Q315r CCD detector. The crystal-to-detector distance was 
279.47 mm, and a total of 360 images were collected with a 
0.5° oscillation angle. The diffraction data set was processed with 
iMOSFLM [37]. In all cases, space group examination was done 
with POINTLESS [39] and intensity scaling and reduction with 
SCALA [39] from the CCP4 suite of programs [40]. Analysis of data 
quality and merohedral twinning was done with both TRUNCATE 
[41] and phenix.xtriage [42]. 

SAD phasing and structure solution 

The structure of LpEstl esterase was determined by single- 
wavelength anomalous diffraction (SAD) at the optimal peak 
wavelength. SAD phasing and model-building was carried out 
with AutoSol and AutoBuild from PHENLX [42] , respectively. AutoSol 
could determine the complete Se substructure and an initial 
electron density map could be calculated. The model automati- 
cally built was made up of 1337 residues in 7 fragments. The 
remaining model was built manually using Coot [43] and the 
refinement was performed with phenix.refine [44]. The data- 
collection and refinement statistics are summarized in Table 1. 
The atomic coordinates and structure factors for LpEstl have 
been deposited in the Protein Data Bank (PDB entry 4c87). 

Structure determination and refinement of the high- 
resolution data set 

Phases for the native and the high-resolution data set were 
obtained by molecular replacement using the program Phaser [45] . 
The atomic coordinates of the SAD structure were used as a 
search model. As above, the model was built manually using Coot 
[40] and the refinement was performed with phenix.refine [44], 
Refinement steps included xyz refinement, TLS, individual atomic 
displacement parameters (ADPs), addition of ligands, and auto- 
matic addition of water molecules using default parameters. The 
complete data-collection and refinement statistics are shown in 
Table 1 . Since a low fraction of twinning was detected in the high- 
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resolution data set (0.17 Britton analysis or H-test; 0.161 
maximum likelihood method), the final refinement was carried 
out by applying the suggested twin law (-h, k, -1). Validation of the 
model was carried out using MolProbity [46]. The atomic 
coordinates and structure factors have been deposited in the 
Protein Data Bank (PDB entries 4c88 and 4c89 for native and 
high-resolution Se-Met labelled LpEstl, respectively). 

Figure preparation 

Model representation was done with PyMOL [47]. 

Analytical chromatography 

Analytical size-exclusion chromatography was performed on a 
Superdex 200 10/300 GL Tricorn column (GE Healthcare) 
equilibrated in 20 mM Tris-HCl, pH 8.0, 0.1 M NaCl, 5 mM 
DTT and 0.04% (w/v) sodium azide. The column was calibrated 
with apoferritin (443 kDa), fS-amylase (200 kDa), alcohol deshy- 
drogenase (150 kDa), bovine serum albumin (66 kDa), ovalbumin 
(45 kDa), carbonic anhydrase (29 kDa), sperm whale myoglobin 
(17 kDa) and vitamin B 12 (1.3 kDa) in the same buffer. The size of 
LpEstl was determined from its K av value (K av = (V e — Vo)/ 
(Vt — Vo); V e : elution volume; V () : void volume; Vp: total volume 
of the column) by interpolation in a calibration semilog plot of the 
molecular mass of the standard proteins versus their K av values. 

Analytical ultracentrifugation 

Equilibrium and velocity ultracentrifugation experiments were 
performed using a Beckman XL-A ultracentrifuge with an An- 
50Ti rotor and standard double sector centrepiece cells. Solvent 
density (1.002 mg ml -1 ) and the partial specific volume of LpEstl 
(0.719) were calculated from the buffer composition (100 mM 
NaCl and 20 mM Tris) and from the predicted amino acid 
composition, respectively, with SEDNTERP [48]. Data from 
sedimentation velocity and equilibrium experiments were analysed 
with the programs Sedfit [49] and Heteroanalysis [50], respec- 
tively. 

Circular dichroism spectroscopy 

Far-UV circular dichroism (CD) measurements were carried out 
on a JascoJ-715 spectropolarimeter equipped with a thermostated 
cell holder and a Peltier temperature control accessory. The 
instrument was calibrated with (+)- 1 0-camphorsulfonic acid. CD 
spectra were recorded in 0. 1 cm path length quartz cells cuvettes 
from 250 to 200 nm at 25 °C, using a protein concentration of 
5.5 uM (1 nm bandwidth, 4 s response, and 20 nm/min scan 
speed). Each spectrum herein presented is the average accumu- 
lation of four scans. Baseline subtraction was performed in all 
cases. Results are expressed as mean residue ellipticity [6]mrw> m 
units of degree cm 2 dmol -1 of amino acid (Mr =110 for this 
protein). Thermal transitions were also analyzed by CD spectros- 
copy by monitoring the variation of the ellipticity at 222 nm as the 
temperature was increased from 20 to 90 °C at 50 °C/h. The 
normalized ellipticity value at each temperature was calculated as 
([6]t - [9]25)/([6]90 - [G]25), where [9] T is the ellipticity value at 
temperature T, and [9] 2 5 and [0] 9O are the ellipticity values at 25 
°C and 90 °C, respectively. Three different samples were analysed, 
although the traces shown correspond to individual samples. 

Enzyme assays and biochemical characterization of 
LpEstl 

Esterase activity was examined spectrophotometrically using p- 
nitrophenyl acetate as substrate. The rate of hydrolysis of pNP- 
acetate for 10 min at 30°C was measured in 50 mM sodium 



phosphate buffer pH 7.0 at 348 nm in a spectrophotometer 
(UVmini-1240 Shimadzu). The amount of protein used was 10 u.g. 
The reaction was stopped by chilling on ice. Controls without 
enzyme were utilized to account for any spontaneous hydrolysis. 
One unit of esterase activity was defined as the amount of enzyme 
required to release 1 umol of ^>-nitrophenol per minute under the 
previously described conditions. 

Substrate specificity of LpEstl was examined using different p- 
nitrophenyl (/>-NP) esters of various chain lengths (C2, C4, C8, 
CI 2, C14 and CI 6) and substrate profile was analysed with a 
library of esters similar to the one reported previously [35]. The 
esters were chosen to identify acyl chain length preferences of 
LpEstl and also its ability to hydrolyse hindered or charged 
substrates. Simple alkyl esters as well as activated esters (vinyl and 
phenyl esters, esters with electron-withdrawing substituents in the 
acyl portion) were included to test whether activated esters would 
react faster. The screening was performed in a 96-well plate Flat 
Bottom (Sarstedt) with a final reaction volume of 200 |il per well, 
each one containing 1 mM substrate in acetonitrile (1% vlti). The 
buffer/indicator solution contained 0.44 mM of />-nitrophenol as 
pH indicator in 1 mM sodium phosphate buffer pH 7.2. Esterase 
(20 ixl of a 0.5 mg ml 1 solution in 1 mM sodium phosphate 
buffer pH 7.2) was added to each well and reactions were followed 
by measuring the decrease in absorbance at 410 nm for 2 h at 
30°C in a Synergy HT BioTek microplate spectrophotometer. 
The incubation time was selected as to maximize the signal-to- 
noise ratio of the absorbance readings. Controls without enzyme 
carried out for each substrate indicated that compounds were 
stable within the time scale of the experiments. Data were 
collected in triplicate and the average activities were quantified. 
Results are shown as means ± standard deviations. 

In order to investigate temperature effect, reactions were 
performed in 50 mM sodium phosphate buffer (pH 7.0) at 4, 20, 
30, 37, 40, 45, 55 and 65 °C. Effect of pH was investigated by 
assaying esterase activity in a range of pH values from 5.5 to 9.0 at 
30 °C. Buffers used were acetic acid-sodium acetate buffer for 
pH 5.5, sodium phosphate buffer for pH 6-7, Tris-HCl buffer for 
pH 8 and glycine-NaOH buffer for pH 9. A 100 mM concentra- 
tion was used in all the buffers. 

For temperature stability measurements, the recombinant 
esterase was incubated in 50 mM sodium phosphate buffer 
pH 7.0 at 20, 30, 37, 45, 55 and 65 °C for 15 min, 30 min, and 
1, 2, 3, 4, 6 and 20 h. After incubation, the residual activity was 
measured as described above. 

Computational methods 

The pdb2pqrpy tool [5 1] was used to estimate the most probable 
protonation states of titratable residues in the protein at a pH of 
7.0 and to add all the missing hydrogen atoms. Hydrogens 
belonging to the catalytic triad (Serl74, His313 and Asp283) were 
manually reoriented in accordance with the known catalytic 
mechanism of the enzyme. 

Affinity potentials [52] within the active site for methyl (sp 3 C), 
hydroxyl (H-bond acceptor/donor sp' 5 O), carbonyl oxygen (H- 
bond acceptor sp 2 O), positively charged amino (H-bond donor 
sp 3 N) and hydrophobic (sp 3 C minus a H-bonding term) probes 
were calculated with our in-house program cGRILL using a grid 
spacing of 0.5 A. Substrates phenyl acetate, triacetin and 
tributyrin were built and assigned point charges and parameters 
from the Merck molecular force field 94 using the openbabel tool. 
Energetically favoured binding poses for these substrates within 
the active site were found by using our in-house CRDOCK tool 
[53]. Briefly, for each substrate a conformational search was 
performed and exhaustive sampling was achieved using the 
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AMBER force field for intermolecular energy evaluation and the 
BFGS algorithm for pose optimization. The final binding modes 
selected were awarded the best score according to the CRScore 
function. 

The nature of the interactions and the residues involved in 
dimerization of LpEstl and its bacterial HSL homologs were 
analyzed by means of the MM-ISMSA method [54], which includes 
the non-bonding term of the AMBER99sb force field [55] and a 
desolvation term as calculated by the Implicit Solvation Method 
(ISM) [56]. The per-residue energy decomposition performed by 
this tool also allowed us to identify those amino acids that contribute 
to the overall stabilization of the dimer. For adaptation to the 
AMBER99sb force field, all proteins were first protonated at 
pH = 7.0 using PDB2PQR, parameterized and geometry optimized 
by following a simple protocol of 500 steps of steepest-descent and 
4500 conjugate-gradient energy minimization with SANDER [57]. 

Supporting Information 

Figure SI Crystals of Se-Met labelled LpEstl. (A) 

Crystallization drop containing both spindle-shaped crystals of 
Se-Met labelled LpEstl, which corresponded to perfectly twinned 
crystals with apparent point group 422 (apparent space group 
7422) and "half crystals resulting from the manipulation of the 
latter, which belonged to the tetragonal 74 space group and did not 
exhibit merohedral twinning. Bar length corresponds to 0.2 mm. 
(B) Diagram explaining the perfect twinning present in the spindle- 
shaped crystals of LpEstl as resulting from the geometrically well 
defined combination of two opposed, untwinned crystals. 
(TIF) 

Figure S2 Tetramers formed by canonical subtype 1 
dimers of enzymes from the hormone-sensitive lipases 
family. (A) tetramer of the thermophilic esterase St-Est from 
Alicyclobacillus acidocaldarius (PDB entry, 3aik). B, tetramer of the 
hyperthermophilic carboxylesterase PestE from the archaea 
Pyrobaculum calidifontis (PDB entry, 2yh2). The dimers at the 
bottom are oriented as in Fig. 3 and are shown as ribbon models, 
whereas the upper dimers are shown as surface plus ribbon models. 
(TIF) 
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